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1.  Problem  Statement 


The  integration  of  real-time  digital  signal  processing  with  sensor  technology  has  spurred  a  renewed 
effort  for  “smart  sensors”  in  military  applications.  For  example,  there  has  been  a  resurgence  in  the  U.S. 
Army’s  interest  in  low  frequency  acoustic  sensors  for  the  identification  and  tracking  of  targets  such  as 
tracked  vehicles,  airborne  vehicles,  and  munition  muzzle  blasts.  For  this  particular  application,  there  exists 
a  need  for  compact,  rugged,  and  very  sensitive  arrays  of  acoustic  sensors.  The  broad  objective  of  this 
research  effort  is  development  of  the  technology  leading  to  a  novel  sensor  based  on  optical  interferometry, 
microelectromechanical  systems  (MEMS),  and  digital  signal  processing  (DSP).  The  synergistic  integration 
of  these  three  technologies  provides  the  advantages  of  high  sensitivity,  frequency  selectively,  and  the 
possibility  of  sensor  array  integration  in  a  compact  and  robust  package.  These  features  are  required  for  use 
on  the  digital  battlefield. 

MEMS  technology  enables  the  fabrication  of  three  dimensional,  miniature  (micron-sized  features),  and 
environmentally  robust  sensing  structures  which  may  be  fabricated  in  a  dense  array  geometry’.  Optical 
interferometry  provides  a  means  of  making  high  sensitivity’  measurements  on  the  MEMS  devices.  The 
measurement  data  immediately  undergoes  digital  signal  processing  in  situ  with  the  sensor  to  extract  the 
desired  information.  Thus,  the  resulting  “smart  sensor”  exhibits  frequency  selectivity  and  high  sensitivity 
in  a  compact  and  robust  package.  These  features  of  the  sensor  enable  digital  signal  processing  techniques 
such  as  frequency  matching  for  the  detection  and  recognition  of  signals  with  distinct  spectral  signatures. 

The  numerous  and  varied  applications  of  this  novel  sensor  include  airborne  acoustic  sensing,  resonant 
magnetic  field  sensing,  ultrasonic  imaging  through  the  ground,  health  monitoring,  meteorological  sensing 
(e.g.,  monitor  and  track  the  progress  of  tornadoes),  non-destructive  evaluation  of  materials,  and  disaster 
relief  (i.e.,  “hearing  through  the  rubble”).  Direct  Army  benefit  from  this  research  effort  will  result  in 
applications  in  the  detection,  classification,  and  triangulation  of  targets  (such  as  artillery,  ground  vehicles, 
and  airborne  vehicles)  from  their  acoustic  signatures.  Another  potential  Army  application  is  magnetic 
sensor  fusing  for  the  targeting  of  massive  metallic  objects  (e.g.,  tanks).  In  addition,  ultrasonic  imaging 
through  the  ground  for  the  detection  of  land  mines  is  yet  another  critical  application  of  the  proposed  smart 
sensor. 

To  achieve  the  desired  smart  sensor  based  on  MEMS  and  optical  sensing,  we  have  identified  two 
specific  technological  areas  of  need:  (1)  development  of  a  robust  interferometric  system  to  measure  the 
out-of-plane  deflection  of  a  MEMS  structure  and  (2)  development  of  DSP  algorithms  for  demodulation  of 
the  interferometer.  The  progress  made  in  each  of  these  two  technological  areas  of  need  is  reported  here. 
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2.  Summary  of  Results 


2.1  Background 

The  rapidly  expanding  field  of  microelectromechanical  systems  (MEMS)  is  experiencing  phenomenal 
growth  in  communication  and  sensing  applications.  MEMS  structures  commercially  available  or  under 
development  for  optical  communication  applications  include  optical  crossconnects  [1],  add/drop 
wavelength  multiplexers  [2],  gain  equalizers  [3],  and  tunable  lasers  [4]  and  filters  [5].  On  the  other  hand, 
MEMS  sensing  applications  have  achieved  commercial  success  in  the  field  microaccelerometers  [6].  Other 
MEMS  sensors  currently  under  investigation  include  pressure  sensors  [7],  magnetic  field  sensors  [8], 
resonant  transducer  sensors  [9],  and  gyroscopic  sensors  [10].  Although  not  all  MEMS  structures  involve 
moveable  components,  the  devices  listed  here  share  the  common  trait  that  they  include  some  out-of-plane 
moveable  component. 

To  develop  viable  MEMS  devices  for  these  sensor  applications,  a  non-invasive  system  for 
characterizing  out-of-plane  displacements  as  a  function  of  the  actuation  parameter(s)  is  desirable.  As  an 
added  benefit,  the  measured  data  frequently  can  be  used  to  deduce  material  property  information. 
Information  of  this  type  provides  feedback  for  adjusting  the  fabrication  and  processing  conditions  to 
achieve  optimal  device  performance.  With  respect  to  device  development,  it  is  also  important  that  the 
sensing  system  be  capable  of  in  situ  characterization.  In  this  context,  it  would  be  desirable  that  the  system 
be  sufficiently  flexible  to  make  measurements  on  MEMS  devices  mounted  on  a  probe  station  or  at  the 
wafer  level  during  process  development  and  manufacturing  [11]. 

The  length  scale  of  out-of-plane  motion  typically  ranges  from  sub-nanometer  to  several  microns  of 
displacement.  The  fine  resolution  and  wide  dynamic  range  afforded  by  optical  detection  techniques 
provides  a  good  match  with  the  demands  of  MEMS  measurements.  As  a  result,  several  optical  systems 
have  been  developed  for  measuring  the  out-of-plane  displacement  of  MEMS  structures.  These  optical 
techniques  may  be  broadly  categorized  as  interferometric  and  non-interferometric.  Non-interferometric 
techniques  typically  modulate  the  power  coupled  into  an  optical  fiber  so  that  it  is  proportional  to  the 
displacement  of  the  MEMS  device.  These  techniques  include  the  optical  beam  deflection  method  [12], 
the  shutter  method,  and  the  lever  method  [13]. 

Interferometric  techniques  can  be  categorized  as  single  point  measurement  techniques  and  full-field 
optical  measurement  techniques.  Single  point  methods  measure  the  transverse  displacement  at  a  single 
point  on  the  MEMS  structure.  Both  bulk  optical  [14]  and  fiber  optic  [15]  interferometers  have  been  used  to 
characterize  the  displacement  at  specific  points  on  the  MEMS  device.  These  approaches  may  be  scanned  to 
yield  linear  or  two-dimensional  displacement  data  [16].  Full-field  optical  measurement  techniques  have 
been  applied  to  MEMS  structures  and  they  include  holographic  interferometry  [17,  18],  moire 
interferometry  [18]  and  stroboscopic  interferometry  [19]. 

Compared  to  sensing  systems  applied  strictly  for  the  development  of  MEMS  devices,  some  sensor 
systems  are  integrated  as  part  of  the  overall  system  architecture.  Detection  mechanisms  used  with  MEMS 
sensors  include  electrostatic,  piezoelectric,  magnetic,  piezoresistive,  and  optical  [20].  Both  integrated 
optical  [21,  22],  and  fiber  optical  [23,  24]  approaches  for  characterizing  MEMS  sensors  have  been  reported 
in  the  literature.  Since  the  fiber  optic  sensing  schemes  are  not  inherently  tied  to  the  MEMS  structure,  they 
are  amenable  to  the  more  general  application  of  in  situ  MEMS  characterization. 

In  this  report,  an  optical  fiber  interferometer  and  digital  demodulation  scheme  are  described  for 
measuring  the  in  situ  deflection  characteristics  of  MEMS  structures.  Common  signal  and  reference  beam 
paths,  together  with  digital  demodulation  techniques,  provide  a  robust  sensing  system,  which  is  tolerant  of 
angular  misalignments  of  the  fiber  probe.  The  system  can  be  expanded  to  multiple  sensor  heads  for  the 
interrogation  of  an  array  of  elements  within  a  MEMS  device.  In  Section  2.2,  the  interferometer  is 
described  with  respect  to  its  optical  configuration  and  digital  demodulation  technique.  The  operational 
characteristics  of  this  system  are  analyzed  in  terms  of  its  dynamic  range  of  measurement  and  its  range  of 
operating  frequencies.  Experimental  results  for  the  characterization  of  a  MEMS  electrostatically 
deflectable  beam  are  also  presented.  In  Section  2.3,  a  novel  DSP-based  demodulation  scheme  is  proposed 
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and  experimentally  verified.  Digital  demodulation  of  the  interferometric  signal  provides  computational 
sophistication,  real-time  operation,  reprogrammability,  and  low  cost.  These  features  are  critical  for  smart 
processing  at  the  sensor  head.  In  Section  2.4  are  the  conclusions. 

2.2  Fiber  Interferometer 

A  fiber-optic-based  interferometer  was  chosen  for  this  application  due  to  its  potential  for  integration 
with  MEMS  structures  and  its  ability  for  distributed  sensing.  In  addition,  the  fiber  interferometer  is 
amenable  to  in  situ  characterization  of  MEMS  structures,  where  MEMS  device  operation  may  be  validated 
at  the  wafer  level  during  fabrication.  To  demonstrate  this  utility  of  the  interferometer,  we  have  integrated 
the  interferometer  with  a  MEMS  probe  station.  In  Section  2.2.1,  we  discuss  the  interferometer’s  optical 
configuration  and  digital  demodulation  scheme.  For  this  system,  a  commercially  available  digital 
demodulator  is  used  to  extract  the  motion  of  the  MEMS  structure.  In  Section  2.2.2,  the  theoretical  range  of 
motion  detectable  by  the  interferometric  system  and  its  frequency  limitations  are  discussed.  Finally  in 
Section  2.2.3,  the  interferometer  is  used  to  characterize  the  motion  of  a  MEMS  flexure  beam.  The  well- 
known  polysilicon  MEMS  cantilever  beam  was  chosen  as  the  test  structure  so  that  mature  modeling 
techniques  could  be  used  to  validate  the  system  response.  The  experimentally  measured  motion  is  in  good 
agreement  with  finite  element  analysis  (FEA)  and  boundary  element  analysis  (BEA)  simulation  of  the 
flexure  beam. 

2.2.1  System  Configuration 

The  interferometer  system  configuration  is  shown  in  Figure  1.  Polarized  light  with  a  wavelength  of 
0.6328  /Jm  is  coupled  into  a  single  mode  optical  fiber.  The  light  is  then  split  with  a  3-dB  fused  fiber 
coupler.  One  of  the  output  pigtails  of  the  coupler  is  attached  to  a  piezoelectric  transducer  stack  (made  from 
lead  zirconate  titanate  or  PZT),  which  moves  the  tip  of  the  fiber  in  a  direction  along  the  axis  of  the  fiber. 
This  pigtail  serves  as  the  fiber  probe  of  the  interferometer.  The  coupler’s  other  pigtail  is  unused.  The 
unused  pigtail  is  immersed  in  index  matching  fluid  to  minimize  back  reflections. 

The  tip  of  the  fiber  probe  is  cleaved  such  that  the  normal  to  the  cleaved  fiber  facet  is  collinear  with  the 
fiber  axis.  The  cleaved  facet  serves  as  both  the  interferometer’s  beam  splitter  and  beam  combiner.  The 
air/ fiber  interface  at  the  cleaved  facet  of  the  fiber  probe  provides  a  4%  power  reflection  back  into  the  single 
mode  fiber.  This  reflected  field  is  the  interferometer’s  reference  beam.  The  light  transmitted  through  the 
fiber  probe’s  facet  reflects  off  the  MEMS  structure  to  be  measured  and  is  coupled  back  into  the  fiber. 
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Figure  1.  Interferometric  System. 
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Figure  2.  Analysis  of  the  power  coupling  efficiency 
of  the  fiber  probe:  (a)  geometry  of  the  system  and 
(b)  power  coupling  efficiency  versus  probe/target 
separation  distance  at  zero  degree  tilt  of  the  fiber  axis 
with  respect  to  the  target  normal. 
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Figure  3.  Analysis  of  the  power  coupling  efficiency 
of  the  fiber  probe:  (a)  geometry  of  the  system  and 
(b)  power  coupling  efficiency  versus  tilt  angle  of  the 
fiber  axis  with  respect  to  the  target  normal  at  a 
separation  distance  of  100  pm. 


Consequently,  the  cleaved  facet  now  acts  as  a  beam  combiner.  The  round  trip  distance  traveled  out  of  the 
fiber,  off  the  MEMS  structure,  and  back  into  the  fiber  forms  the  signal  path  of  the  interferometer.  The 
reference  and  signal  fields  co-propagate  back  through  the  fiber  until  they  reach  the  fused  fiber  coupler, 
where  an  additional  3-dB  power  reduction  is  encountered.  Finally,  the  interference  signal  is  detected  with  a 
photodetector. 

The  coupling  efficiency  of  the  signal  field  back  into  the  fiber  probe  is  a  critical  feature  of  this 
interferometer  configuration.  We  have  modeled  the  coupling  efficiency  of  the  fiber  probe,  assuming 
Gaussian  field  propagation  out  of  the  fiber  and  reflection  off  a  tilted  specular  surface.  The  power  coupling 
efficiency  versus  the  probe-to-target  separation  distance  is  plotted  in  Figure  2.  Given  zero  tilt  of  the  fiber 
probe  with  respect  to  the  target  surface  normal,  the  probe-to-target  separation  distance  must  be  less  than 
195  /im  for  greater  than  4%  power  coupling  of  the  signal  field  back  into  the  fiber. 

Figure  3  shows  the  power  coupling  efficiency  versus  probe  tilt  angle  at  a  separation  distance  of 
100  //m  .  From  an  implementation  perspective,  high  tolerance  to  angular  misalignment  of  the  probe  is 

desired.  The  bare  fiber  probe  tolerates  approximately  ±  4°  of  misalignment  at  a  separation  distance  of 
100  //in  .  This  limitation  is  sufficient  for  the  in  situ  characterization  of  MEMS  structures. 

The  optical  configuration  described  above  has  been  used  with  varying  demodulation  schemes  to 
characterize  MEMS  sensors  [15],  hard  disk  surfaces  [25],  and  biological  membranes  [26].  The  common 
signal  and  reference  paths  in  this  optical  configuration  eliminate  signal  fading  due  to  non-signal  induced 
polarization  drift  and  optical  path  length  drift  in  the  fiber.  Also,  the  common  signal  and  reference  paths 
obviate  the  need  for  high-cost  polarization-preserving  single  mode  fiber. 
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Figure  4.  Interferometric  measurement  of  a  PZT-driven  mirror.  The 
upper  oscilloscope  trace  is  the  signal  driving  the  PZT-mounted  mirror  and 
the  lower  trace  is  the  amplified  output  signal  from  the  photodetector.  The 
output  signal  exhibits  the  characteristic  waveform  for  a  4 K  phase  shift. 


Demodulation  of  the  optical  signal  is  implemented  by  sampling  the  photodetected  signal  and 
performing  digital  signal  processing  on  it.  The  OPD-200  Digital  Demodulator  is  a  commercially  available 
instrument  from  Optiphase™,  Inc.,  which  was  used  to  perform  the  digital  demodulation  [27].  The 
demodulation  algorithm  is  based  on  a  passive  homodyne  scheme  using  a  phase  generated  carrier  (PGC) 
[28].  In  our  interferometer,  the  PGC  is  imposed  on  the  signal  by  modulating  the  open  air  signal  path 
between  the  fiber  probe  and  the  MEMS  structure. 

We  verified  the  digital  demodulation  of  the  fiber  interferometer  by  comparing  the  phase  output  from 
the  OPD-200  with  a  direct  measurement  of  the  phase.  A  PZT-mounted  mirror  was  used  in  place  of  the 
MEMS  structure  in  Figure  1.  The  mirror  was  sinusoidally  vibrated  at  1  kHz  with  varying  drive  voltages 
applied  to  its  PZT  stack.  The  direct  measurement  of  the  mirror  motion  was  obtained  by  noting  that  at 
integer  multiples  of  In  radians  in  the  interferometer’s  optical  path  difference,  characteristic  waveforms  are 
produced  by  the  interfering  beams.  For  example,  the  bottom  trace  in  Figure  4  illustrates  the  photodetected 
waveform  for  a  An  radian  optical  path  length  difference.  Direct  measurement  data  corresponding  to 
optical  path  length  differences  of  2 n  and  An  versus  the  voltage  applied  to  the  PZT  stack  are  plotted  as 
X ’s  in  Figure  5.  Agreement  between  the  direct  measurement  result  and  the  demodulation  output  from  the 
OPD-200  (plotted  as  o’s)  is  evident  in  Figure  5. 


Figure  5.  Comparison  of  the  OPD-200  digital 
demodulator  output  with  a  direct  measurement 
of  the  interferometer  phase:  “  X  ”  =  direct 
measurement  data  point;  “On  =  OPD-200  data 
point;  and  solid  curve  =  linear  fit  to  the  OPD- 
200  data. 
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2.2.2  System  Dynamic  Range 

The  displacement  measurement  and  frequency  response  of  the  interferometer  system  described  above 
has  a  large  dynamic  range.  The  digital  demodulation  technique  facilitates  both  fractional  fringe  and  fringe 
counting  interrogation  of  the  interferometric  signal.  This  technique  has  a  theoretical  dynamic  range  greater 

8 

than  10  for  measuring  surface  displacements. 

The  ultimate  system  resolution  of  a  displacement  measuring  interferometer  refers  to  a  limit  in  the 
system’s  ability  to  detect  surface  changes.  This  limit  is  expressed  in  terms  of  the  smallest  displacement 
that  can  be  detected  if  the  signal-to-noise  ratio  (SNR)  is  unity  [29].  Photodetector  quantum  noise  and 
analog-to-digital  converter  quantization  noise  are  examples  of  noise  sources  that  cannot  be  eliminated  from 
the  system.  In  the  present  system,  the  detector  quantum  noise  (also  called  shot  noise)  sets  the  limit  to  the 
ultimate  performance  of  the  interferometric  system.  We  have  determined  the  quantum-noise-limited 
sensitivity  of  the  interferometer  shown  in  Figure  1  to  be  24  pm.  The  derivation  of  the  ultimate  system 
resolution  for  the  classical  Michelson  interferometer  is  given  in  Appendix  A. 

The  degree  to  which  this  ultimate  resolution  may  be  achieved  depends  on  numerous  factors.  The 
foremost  consideration  for  this  interferometer  configuration  is  the  field  contrast  ratio.  In  practice,  it  is  very 
difficult  to  secure  a  unity  field  contrast  ratio.  For  values  of  field  contrast  not  equal  to  one,  the  fringe 
contrast  and  the  SNR  decrease,  thereby  increasing  the  minimum  detectable  displacement.  In  practice,  we 
have  achieved  a  system  resolution  of  several  nanometers. 

The  number  of  bits  in  the  digital  demodulator,  which  are  used  to  keep  the  fringe  count,  limits  the  full 
measurement  range  of  the  interferometer.  This  limitation  places  an  upper  bound  of  2.6  mm  on  the 
measurement  range.  The  practical  measurement  range  of  several  nanometers  to  more  than  two  and  a  half 
millimeters  is  well-suited  to  the  out-of-plane  characterization  of  MEMS  structures. 

The  interferometric  system  is  capable  of  measuring  displacements  ranging  in  frequency  from  dc  to 
one-half  the  frequency  of  the  phase  generated  carrier.  The  OPD-200  provides  the  drive  signal  for  the  phase 
generated  carrier  and  is  capable  of  modulation  rates  up  to  95  kHz.  In  our  interferometric  configuration, 
however,  the  PZT  stack  limits  the  modulation  frequency.  The  PZT  and  drive  electronics  in  our  system 
limit  the  carrier  frequency  to  10  kHz.  Therefore,  our  system  is  capable  of  measuring  displacements  ranging 
from  dc  to  5  kHz. 

2.2.3  Characterization  of  a  MEMS  Flexure  Beam 

A  MEMS  flexure  beam  is  a  clamped-free  cantilever  beam,  where  a  small  section  of  the  beam  material 
near  the  base  of  the  post  has  been  removed  to  produce  a  hinge  point  on  the  beam  (see  Figure  6). 
Polysilicon  flexure  beams  with  dimensions  given  in  Figure  6  were  fabricated  via  the  Multi-User  MEMS 
Processes  (MUMPs®)  foundry  service  available  from  JDS  Uniphase  [30].  We  used  the  fiber 


Figure  6.  Geometry  of  a 
MEMS  flexure  beam:  (a) 
top  view  and  (b)  side  view. 
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Figure  7.  Interferometric  measurement  of  a  MEMS  flexure  beam.  The 
upper  oscilloscope  trace  is  the  signal  driving  the  MEMS  structure  and  the 
lower  trace  is  the  OPD-200  demodulator’s  output.  The  demodulated 
output  signal  scale  corresponds  to  25.18  nm/V. 


interferometer  to  measure  the  transverse  deflection  of  the  flexure  beam  under  electrostatic  actuation. 
Additionally,  we  simulated  the  motion  of  the  flexure  beam  using  a  coupled  BEA  and  FEA  software  tool 
(IntelliSuite™  Windows  NT  Version  5.1)  available  from  Corning  IntelliSense.  In  the  coupled  approach, 
BEA  uses  a  surface  mesh  and  solves  the  electrostatic  aspect  of  the  problem,  while  FEA  uses  a  volume  mesh 
and  solves  the  mechanical  aspect  of  the  problem.  When  considering  electrostatically  actuated  devices,  the 
coupled  approach  tends  to  several  advantages  over  a  pure  finite  element  approach  including  faster 
simulation  times  and  higher  accuracy  results. 

The  flexure  beam  was  electrostatically  excited  into  transverse  motion  with  the  application  of  a  2.5  Vo-pk 
sinusoid  oscillating  at  1  kHz.  An  example  of  the  interferometer’s  measurement  is  shown  in  Figure  7. 
Channel  1  is  the  upper  trace  and  shows  the  voltage  signal  driving  the  flexure  beam.  Channel  2  is  the  lower 
trace  and  shows  the  demodulation  output  from  the  OPD-200.  Given  the  round  trip  path  taken  by  the  signal 
beam  in  our  interferometer  configuration  (see  Figure  1),  the  peak-to-peak  deflection  of  the  device  under 
test  may  be  described  as 


Mpk-pk 


A-A0 

4n 


where  A(Z>  is  the  phase  difference  between  the  signal  and  reference  beams.  The  digital  demodulator  scale 
is  set  to  0.5  radian /V,  which  corresponds  to  a  peak-to-peak  beam  deflection  of  5.5  nm. 

Figure  8  shows  the  numerical  simulation  of  the  flexure  beam’s  deflection  based  upon  the  coupled  BEA 
and  FEA  technique.  The  mesh  size  was  30  pm  and  the  following  material  properties  were  used  for  the 
polycrystalline  silicon  beam: 

Young’s  Modulus:  160  GPa 

Density:  2.3  g/cm3 

Poisson’s  Ratio:  0.226 
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Displacement  Z  vs  Time 


IntelliSuite  PC  Version,  IntelliSense  Corporation (c) 


Figure  8.  Coupled  boundary  element  analysis  and  finite 
element  analysis  simulation  of  a  MEMS  flexure  beam.  The 
signal  driving  the  MEMS  beam  is  2.5sin(2;r  ■  lOOOt)  V. 


The  data  in  Figure  8  indicates  that  the  flexure  beam  deflects  toward  the  substrate  with  a  frequency  of  2  kHz 
and  a  peak-to-peak  amplitude  of  4.0  nm.  This  simulation  result  compares  favorably  with  the  experimental 
data. 

Of  particular  interest  is  the  oscillation  frequency  of  the  flexure  beam.  The  demodulation  output  shown 
in  the  lower  trace  of  Figure  7  indicates  that  the  flexure  beam  is  oscillating  at  2  kHz  or  twice  the  excitation 
frequency.  During  the  negative  voltage  portion  of  the  excitation  sinusoid’s  period  (see  Figure  9a),  the  E 
field  is  directed  in  the  (+z)  -direction.  The  negative  charge,  —Q  ,  accumulated  on  the  beam’s  electrode 

experiences  a  force  given  by  F  =  —QE ;  so,  the  beam  is  deflected  towards  the  substrate.  During  the 

positive  voltage  portion  of  the  excitation  sinusoid’s  period  (see  Figure  9b),  the  E  field  is  directed  in  the 
(— z)  -direction.  The  positive  charge,  Q ,  accumulated  on  the  beam’s  electrode  experiences  a  force 

F  =  QE  ;  so,  the  beam  is  again  deflected  towards  the  substrate. 

Electrostatically  driven  MEMS  structures  experience  only  electrostatic  attraction  (as  opposed  to 
electrostatic  repulsion).  Thus,  for  a  bipolar  excitation  signal  with  zero  mean,  the  oscillation  frequency  of 
the  beam  is  twice  that  of  the  excitation  signal.  This  frequency  doubling  effect  may  be  removed  by  adding  a 
dc  bias  to  the  excitation  signal  such  that  the  excitation  signal  is  always  a  positive  voltage  or  always  a 
negative  voltage.  Figure  10  shows  a  biased  excitation  signal  (upper  trace)  and  the  raw  interferometer 
output  signal  prior  to  digital  demodulation  (lower  trace).  The  excitation  signal  is  5  + 5sin(2;r  •  1000?)  V. 
Note  that  both  the  excitation  and  the  raw  interferometer’s  response  are  oscillating  at  1  kHz. 

Two  points  of  clarification  are  necessary  with  respect  to  the  interferometer  data  displayed  in  Figure  7. 
First,  the  interferometer  measures  relative  phase  shifts  between  a  reference  light  wave  and  a  signal  light 
wave.  The  demodulator  reads  both  the  offset  phase  (comprising  the  optical  path  mismatch)  as  well  as  the 
dynamic  signal  created  by  the  beam  flexure.  Therefore,  a  zero  voltage  reading  at  the  demodulator  output 


zu 


Figure  9.  Electrostatic  attraction  of  a  MEMS  cantilever 
beam:  (a)  negative  voltage  applied  to  the  beam  and  (b) 
positive  voltage  applied  to  the  beam. 


Figure  10.  Removal  of  the  frequency  doubling  effect.  The  upper  oscilloscope 
trace  is  a  biased  excitation  signal  driving  the  MEMS  structure  and  the  lower 
trace  is  the  raw  interferometer  output  signal  prior  to  demodulation.  The  signal 
driving  the  MEMS  beam  is  5  +  5  sin(2;r  •  1  OOOt)  V. 
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does  not  necessarily  correspond  to  zero  deflection.  Second,  an  artifact  of  the  digital  demodulation  scheme 
is  a  fixed  time  delay  in  the  output  data,  which  corresponds  to  the  time  for  one  demodulation  cycle.  For  the 
modulation  frequency  used  in  Figure  7,  the  delay  is  about  124  jUs  . 

Finally,  we  measured  the  deflection  at  the  tip  of  the  flexure  beam  as  a  function  of  the  excitation 
voltage  at  1  kFlz.  As  shown  in  Figure  11,  the  experimentally  measured  data  is  in  reasonably  good 
agreement  with  the  simulation  data.  Flowever,  a  discrepancy  in  deflection  amplitudes  at  the  larger 
excitation  voltages  exists  between  the  experimental  and  the  simulation  data.  Further  simulation  results 
indicate  that  magnitude  of  this  discrepancy  is  not  attributable  to  practical  variations  in  the  beam  parameters 
(e.g.,  Young’s  modulus  and  beam  length)  or  to  the  onset  of  spontaneous  collapse  of  the  flexure  beam. 
Furthermore,  simulation  indicates  that  the  flexure  beam  has  a  transverse  mode  resonance  at  73  kFlz.  Our 
beam  oscillation  frequency  of  2  kFlz  is  sufficiently  removed  from  resonance  such  that  overshoot  of  the 
beam  cannot  be  the  cause  of  the  discrepancy.  This  discrepancy  between  the  experiment  and  simulation  is 
currently  under  investigation. 


Figure  11.  Tip  deflection  of  the  MEMS  flexure  beam  versus 
applied  voltage  at  a  frequency  of  1  kFlz:  “  X  ”  =  FEA 
simulation  data  point;  solid  curve  =  quadratic  curve  fit  to  the 
simulation  data;  “0”  =  experimental  data  point;  and  dotted 
curve  =  quadratic  curve  fit  to  the  experimental  data. 
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2.3  A  Novel  Digital  Demodulation  Algorithm 

An  optical  interferometer  may  be  utilized  to  measure  the  displacement  of  an  object  from  a  reference 
location  by  observing  the  interference  pattern  between  two  beams  of  light  that  vary  in  phase  as  caused  by 
the  displacement  of  the  object.  Demodulation  of  the  interferometer  refers  to  extracting  the  motion  of  the 
object  from  information  contained  within  the  optical  interference  pattern.  In  this  section,  we  describe  a 
novel  interferometric  demodulation  algorithm  implemented  with  a  digital  signal  processor.  Digital  signal 
processing  (DSP)  is  a  useful  method  of  analyzing  and  manipulating  an  interferometer  signal.  The 
architecture  and  instruction  set  of  digital  signal  processors  are  optimal  for  the  high-speed,  mathematically 
intensive  calculations  involved  in  demodulation.  The  advantages  of  DSP  include  mathematical 
sophistication,  reprogrammability,  and  cost-effectiveness.  In  our  demodulation  algorithm,  we  assume  that 
the  target  is  confined  to  sinusoidal  motion  of  a  known  frequency.  This  algorithm  therefore  may  be  applied 
to  systems  vibrating  at  a  known,  fixed  frequency,  including  transducer  calibration  [14],  acoustic  sensing 
[31],  and  microelectromechanical  systems  (MEMS)  characterization  [32].  Our  algorithm  is  based  on  use 
of  a  PGC  [28].  The  frequency  components  of  the  resultant  signal  from  the  PGC  and  the  vibrating  target 
contain  information  concerning  the  mean  optical  path  length  difference,  the  modulation  depth  of  the  carrier, 
and  the  vibration  amplitude  of  the  target.  These  frequency  components  may  be  manipulated  to  extract  only 
the  desired  target  vibration  amplitude.  In  Section  2.3.1,  the  algorithm  and  its  theoretical  dynamic  range  of 
measurement  are  discussed.  In  Section  2.3.2,  experimental  results  are  presented  to  verify  the  algorithm. 

2.3.1.  Dynamic  Range  of  Measurement 

The  application  of  this  interferometric  demodulation  algorithm  is  the  measurement  of  the  vibration 
amplitude  of  a  structure  given  its  frequency.  Figure  12  shows  a  Michelson  interferometer  equipped  with  a 
reference  mirror  that  is  modulated  at  a  known  frequency  cor  and  an  amplitude  b  =  A/4,  where  A  is  the 
wavelength  of  the  laser  light.  A  target  is  assumed  to  be  vibrating  sinusoidally  at  a  known  frequency  to,  that 
is  much  greater  than  cor.  The  photodetector  signal  contains  information  concerning  the  target  vibration 
amplitude  a.  A  DSP  processor  has  been  incorporated  to  provide  the  feedback  necessary  for  stabilization,  as 
well  as  to  demodulate  the  interferometer  output  signal  and  to  extract  the  target  vibration  amplitude  a. 

The  ac  component  of  the  output  optical  power  of  the  interferometer  is  given  by 

Pout,ac  K  cos[d>  E  +  2ak  sin  (ft;,  f  +  </>,)+ 2bk  Am(co,t  +  <j)r )],  (1) 

where  K  is  the  system  gain  constant;  <D£  is  the  equilibrium  optical  path  length  difference;  k  is  the 
propagation  constant  of  light  in  air;  and  <f>r  and  (f),  are  the  phases  of  the  reference  and  target  beams, 


Figure  12.  Michelson  interferometer 
with  feedback  controller. 
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respectively.  As  the  target  vibration  amplitude  increases,  the  interferometer  output  signal  exhibits  higher- 
order  harmonics.  By  applying  the  Fourier-Bessel  expansion  [33]  to  Equation  (1),  the  harmonic  content  of 
this  signal  is  revealed: 

Pouuac  “  K{co^E)j0{2ak)J0{2bk) 

-  2sin(<J>£)j0(2a£),/1  {2bk)s'm(cort  +  <j>r ) 

+  2  cos(d>£  )J0  (2  ak)j2  (2bk)s'm(2cort  +  2<f>r ) 

-2sin(<&E)jo{2ak)j}{2bk)sin(3(0rt  +  3^,.) 

-  2  sin(<J>£  )j,  (2  ak)j0  {2bk)s,'m{colt  +  <f>, ) 

-  2  cos(<J>£  )j,  ( 2ak)j]  (2bk)cos(colt  +  cort  +  </>,+  </ ) 

-  2  cos(<J>£ )j, (2 ak)j]  (2bk)cos(-  co,t  +  cort  -</>,+  </>,. ) 

+  2cos(4>£)y2  (2aA:)J0(2M)sin(2ft;/t  +  2(Z>() 

-  2  sin(<J>£  )j2  (2ak)jl  (2M)  sin(2ry,  t  +  cort  +  2q \  +  (f. ) 

-  2  sin(<J>£ )j2 {2ak)jx (2M)  sin(-  2 a>,t  +  cort  -  2(f),  +  <f>r).. .}, 

where  /„(•)  is  an  nth  order  Bessel  function  of  the  first  kind.  Note  that  many  of  the  harmonics  contain 
identical  factors  in  their  amplitudes.  Two  carefully  selected  frequency  components  may  be  detected  and 
manipulated  for  demodulation.  These  signal  components  are  expressed  as 

Xc (a)=  K* -2 sin(<I>£ )j0 (2 ak)jt ( 2bk)sin(a>rt  +  (j)r )  and 

XD  (a)  =  K*  -2  sin(<f>£  )j2  (2ak)jt  (2  bk)  sin(2  cott  +  cort  +  2  (f>t  +</>,.)■  (3) 

The  angular  frequencies  of  Xc  and  XD  are  cor  and  2 co,  +  cor,  respectively.  These  specific  terms  have 
been  chosen  such  that  sin(<f>£),  J\{2bK),  and  K  are  strategically  removed  by  taking  the  ratio  of  their 
amplitudes.  This  reveals  a  measurable  quantity  that  is  only  a  function  of  the  desired  target  vibration 
amplitude  a,  given  by  the  equation 


|XD(a)|  |T2(2aA:)| 

This  is  significant  in  that  it  allows  demodulation  to  be  completely  independent  of  the  equilibrium  optical 
path  length  difference  <f>£,  the  modulation  depth  2b k.  and  the  system  gain  K.  By  implementing  a  root- 
finding  algorithm,  such  as  bracketing  and  bisection  [34],  the  argument  2 ak  of  Equation  (4)  may  be 
extracted.  Further  division  by  2k  reveals  the  target  vibration  amplitude,  a. 

Figure  13  plots  Equation  (4)  versus  amplitude  a,  indicating  a  maximum  detectable  amplitude.  Because 
the  demodulation  scheme  involves  measuring  the  value  of  a  function  that  has  reoccurring  poles  and  zeros 
with  increasing  values  of  a,  amplitude  certainty  can  only  be  guaranteed  below  the  first  zero  of 
\Ja{2ak)\l\J2{2ak)\.  This  is  known  as  the  distance  ambiguity  function.  For  an  operating  wavelength  X  = 
632.8  nm,  amax  ~  120  nm.  The  minimum  detectable  displacement  due  to  quantum  shot  noise  for  this  phase¬ 
generated  carrier  demodulation  scheme  is  amm  =  5  pm  (see  Appendix  A).  Therefore,  the  measurement 
dynamic  range  of  this  system  is  on  the  order  of  105. 

2.3.2.  Experimental  Verification 

The  TMS320C31  32-bit  floating-point  digital  signal  processor,  provided  with  the  Digital  Signal 
Processor  Starter’s  Kit  (DSK)  available  from  Texas  Instillments,  Inc.,  was  used  to  stabilize  the 
interferometer’s  operating  point  at  quadrature  and  to  demodulate  the  interferometer  output  signal.  The 
TMS320C31  DSK  was  programmed  using  the  digital  stabilization  scheme  based  on  the  Goertzel  algorithm 
(GA)  developed  by  our  group  [35].  The  Goertzel  algorithm  is  a  filtering  adaptation  of  the  discrete  Fourier 
transform  (DFT).  Compared  to  the  fast  Fourier  transform  (FFT),  the  Goertzel  algorithm  is 
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Distance  Ambiguity  Function,  A  =  632.8  nm 


Vibration  Amplitude,  a  [nm] 

Figure  13.  The  distance  ambiguity  function,  |^/q(2«A:)| /|*/2 )| ,  for 
a  fractional  fringe  interferometer.  The  wavelength  A  detennines  the 
maximum  amplitude  flmax  that  can  be  unambiguously  demodulated. 


computationally  efficient  when  a  low  number  of  frequencies  are  to  be  detected.  An  algorithm  was  then 
developed  to  perform  the  demodulation  of  an  interferometer  signal  and  the  extraction  of  the  target 
amplitude  proposed  in  the  previous  section.  The  structure  of  the  program  is  outlined  in  Figure  14. 


Figure  14.  Flow  chart  for  the 
demodulation  of  an  optical 
interferometer  using  the  Goertzel 
algorithm  (GA). 
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A  test  of  the  demodulation  process  involved  creating  a  signal  that  sufficiently  represented  the 
interferometer  output  by  simulating  only  those  frequency  components  necessary  for  demodulation.  Since 
Equation  (2)  indicates  that  the  interferometer  output  signal  is  composed  of  a  sum  of  harmonics,  a 
simulation  of  the  interferometer  output  signal  was  synthesized  by  summing  a  pair  of  sinusoids.  Recalling 
Equation  (4),  the  simulated  target  vibration  amplitude  a  was  varied  by  adjusting  the  ratio  of  the  sinusoid 
amplitudes. 

The  demodulation  parameters  were  set  for  a  reference  frequency  of fr  =  315  Hz  and  a  target  frequency 
ft  =  2008  Hz.  Therefore,  the  frequencies  to  be  detected  were  fc  =  315  Hz  and  fD  =  4331  Hz.  A 
multifunction  synthesizer  was  used  to  create  the  simulated  interferometer  signal 

v(t)  =  \XC | sin(2;r  ■  3 1 5  •  t)  +  \XD | sin(2;r  • 43  3 1 • 7 )  [V].  (5) 

A  set  of  simulated  signals  was  defined  according  to  the  parameters  listed  in  Table  1.  The  range  of  the 
analog-to-digital  converter  (ADC)  of  the  DSK  was  limited  to  ±1.5  V.  The  magnitudes  of  Xc  and  XD  were 
selected  such  that  their  sum  did  not  exceed  1.5  V.  Table  2  and  Figure  15  show  that  the  DSK  successfully 
demodulated  the  simulated  interferometer  output  signal  precisely  to  within  0.1  nm  and  accurately  to  within 
1%  of  the  simulated  value. 


The  next  objective  was  to  test  the  stabilization/demodulation  routine  on  an  appropriate  system  that  is 
easily  characterized.  To  accomplish  this  task,  a  PZT  stack  was  chosen  to  be  the  vibrating  target.  The 
optical  interferometer  configuration  is  described  in  Section  2.2.1  and  shown  in  Figure  1. 

The  distance  through  which  a  vibrating  object  oscillates  may  be  estimated  by  observing  the 
interference  pattern  created  in  the  interferometer  output  signal.  Peak-to-peak  phase  shifts  at  multiples  of  2n 
display  characteristic  waveforms  that  are  easily  identifiable  when  viewed  with  an  oscilloscope.  When 
operating  at  quadrature,  these  waveforms  are  symmetric  about  their  average  as  illustrated  in  Figure  16. 
According  to  this  phenomenon,  the  PZT  was  characterized  by  a  direct  measurement  of  phase.  The  voltages 
corresponding  to  the  target  PZT  zero-to-peak  phase  shifts  2 ak  =  \  K  2k,  3 k}  were  recorded.  The  phase 
shifts  were  converted  to  amplitudes  through  the  relationship 


a  = 


(6) 


Simulated 

\Xc{a 

)|  \la{2uk\ 

\XC\  [mV] 

\XD\  [MV] 

Amplitude,  a  [nm] 

)|  Vl^ak} 

5 

810.11926 

1000 

1.234 

10 

201.52956 

1000 

4.962 

20 

49.38136 

1000 

20.25 

30 

21.20452 

1000 

47.16 

40 

11.34114 

1000 

88.17 

50 

6.77411 

1000 

147.6 

60 

4.29140 

1000 

233.0 

70 

2.79238 

1000 

358.1 

80 

1.81725 

500 

275.1 

90 

1.14629 

500 

436.2 

100 

0.66371 

500 

753.3 

110 

0.30370 

303.7 

1000 

120 

0.02657 

26.57 

1000 

Table  1 .  The  values  \XC\  and  \XD\  used  to  simulate  an  interferometer 
signal  corresponding  to  various  target  amplitudes  a. 
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Measured  Vibration  Amplitude,  a  [nm] 


Simulated 

Average  Measured 

Probable  Error  for  24 

%  Error 

Amplitude,  a  [nm] 

Amplitude,  a  [nm] 

samples  [nm] 

5 

4.952 

8.972E-02 

-0.9573% 

10 

9.909 

4.360E-02 

-0.9097% 

20 

19.97 

2.533E-02 

-0.1514% 

30 

29.94 

1.994E-02 

-0.2022% 

40 

39.95 

1.628E-02 

-0.1174% 

50 

49.95 

0 

-0.0901% 

60 

59.95 

0 

-0.0914% 

70 

69.94 

0 

-0.0924% 

80 

79.98 

0 

-0.0197% 

90 

89.97 

0 

-0.0284% 

100 

100.0 

1.941E-02 

0.0028% 

110 

110.0 

0 

0.0124% 

120 

120.0 

0 

0.0031% 

Table  2:  Simulation  results.  The  DSK  successfully 

computed  the  simulated  amplitude  with  less  than  1%  error. 


C31  DSK  Demodulation:  Measured  vs.  Simulated 
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Figure  15.  Plot  of  the  simulation  results. 
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Characteristic  Phase  Waveforms 


2  ak  =  tt 


Time 


2ak  =  2tt 


Figure  16.  Characteristic  waveforms  at  the  output  of  an  interferometer 
stabilized  at  quadrature.  The  reference  arm  is  static  and  the  target  is  vibrating 
at  an  amplitude  2 ak. 


These  known  amplitudes  are  plotted  versus  their  respective  driving  voltages  in  Figure  17.  A  linear  fit 
to  this  data  reveals  a  non-zero  a  when  0  V  is  applied.  Therefore,  the  voltage-amplitude  relationship  for  the 
PZT  was  created  in  a  piece-wise  linear  fashion  such  that 

,  ^(2at{v)YX\_  f  0.7812F,  2 ak<n  (?) 

°  In  A4j~{o.8370F-12.1094,  2ak  >  n 

where  a  is  in  units  of  [nm]  and  V is  in  units  of  [mV]. 


Figure  17.  Linear  approximation 
of  the  V0_pk  —  a  characteristic  of 
a  vibrating  PZT  stack. 
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\J0(2ak)\ 

\J2{2ak)\ 


Figure  18.  Plot  of  the  theoretical  \Jo(2ak)  /  and  the 

experimental  |X^(a)|/|X£)(a)|  as  a  function  of  a. 


Equation  (7)  was  then  used  to  test  the  ability  of  the  DSK  to  compute  the  demodulation  amplitude  a 
from  the  two  signals  \Xc{a)\  and  \XD(a)\.  The  code  was  written  to  return  24  values  of  \Xc(d)\l\XD{d)\.  Using 
Equation  (7)  and  varying  the  driving  voltage  V.  the  vibration  amplitude  of  the  PZT  was  varied  between  0 
and  450  nm.  Figure  18  illustrates  a  plot  of  experimental  \Xc(a)\l\XD{a)\  versus  a  superimposed  with 
theoretical  \J0(2ak)\/\J2(2ak)\  versus  a.  Figure  18  shows  that  the  measurement  of  \J0(2ak)\/\J2(2ak)\  is 
accurate  over  a  large  range  of  amplitudes  provided  that  the  corresponding  amplitude  a  is  not  located  near  a 
pole  or  a  zero,  where  small  deviations  in  a  cause  large  changes  in  the  value  of\J0{2ak)\l\J2{2ak)\. 

The  demodulation  algorithm  assumes  that  a  target  is  vibrating  at  an  amplitude  below  the  first  zero  of 
the  distance  ambiguity  function.  Figure  18  suggests  that  if  the  target  vibration  amplitude  is  known  to  be 
located  between  consecutive  poles  and  zeroes,  the  first  zero  no  longer  remains  a  limit.  Currently,  however, 
the  distance  ambiguity  function  is  a  major  limitation  of  the  demodulation  algorithm.  The  first  zero  of 
\J0{2ak)\l\J2{2ak)\  occurs  at  2 ak  ~  2.4,  or  for  A  =  632.8  nm,  amax  =120  nm. 

The  demodulation  parameters  used  for  simulation  were  recalled  for  this  test.  For  a  reference  frequency 
of  fr  =  315  Flz  and  a  target  frequency  f  =  2008  Flz,  the  frequencies  to  be  detected  were  /c  =  315  Flz  and  fD  = 
433 1  Flz  as  specified  by  Equation  (3).  The  demodulation  program  was  executed,  and  the  results  are  plotted 
in  Figure  19  versus  the  amplitude  predicted  by  Equation  (7).  Figure  19  indicates  that  for  amplitudes  a  >  60 
nm,  the  data  obtained  from  DSP  demodulation  closely  matches  the  analytical  data.  For  a  <  60  nm,  the 
DSK  measures  approximately  5  nm  below  the  predicted  value  of  a.  These  deviations  are  attributed  to  non¬ 
linear  operation  of  the  PZT  at  low  drive  voltages. 

To  verify  the  accuracy  of  the  measurement  and  to  confirm  the  non-linearity  of  the  PZT  stack,  the  OPD- 
200™  Digital  Demodulator  by  Optiphase,  Inc.  [25]  was  implemented  in  place  of  the  DSK.  The  OPD- 
200™  is  a  commercially  available  time-domain  demodulation  tool  useful  in  interferometric  measurement. 
The  demodulation  algorithm  is  based  on  a  passive  homodyne  scheme  using  a  PGC  modulation.  The  OPD- 
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DSP-Demodulated  Amplitude  vs.  Predicted  Amplitude 
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Figure  19.  Plot  of  the  amplitude  measured  using  the  DSP 
technique  versus  the  amplitude  predicted  by  Equation  (7). 


200™  provides  the  drive  signal  for  the  PGC  and  is  capable  of  modulation  rates  lip  to  95  kHz.  The  OPD- 
200™  can  measure  arbitrary  displacements  as  well  as  vibrations  ranging  in  frequency  from  dc  to  one -half 
the  frequency  of  the  phase-generated  carrier. 

The  available  PZT  and  drive  electronics  in  our  system  limited  the  carrier  frequency  to  3.5  kHz.  The 
OPD-200™  was  used  to  demodulate  the  interferometer  signal  with  a  target  PZT  oscillating  at  512  Hz.  Our 
DSP-based  demodulation  system  was  used  to  demodulate  the  interferometer  signal  with  the  target  PZT 
oscillating  at  1024  Hz.  This  discrepancy  between  target  frequencies  of  the  two  systems  was  assumed  to  be 
acceptable.  Figure  20  illustrates  the  comparison  of  the  measurements  made  with  the  DSP  demodulation 
scheme  versus  those  obtained  using  the  OPD-200™.  The  experimental  data  are  nearly  identical,  thereby 
verifying  our  DSP  demodulation  system  over  the  desired  measurement  range. 


Verification  of  the  DSP  Demodulation  Using  the  OPD-200 
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Figure  20.  Plot  of  the  amplitude  measured  using  the  DSP 
technique  versus  the  amplitude  measured  using  the  OPD- 
200™  Digital  Demodulator. 
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2.4  Conclusion 


We  have  reported  a  fiber  optic  interferometer,  which  is  suitable  for  the  in  situ  characterization  of 
MEMS  structures.  The  optical  configuration  of  the  interferometer  is  robust  with  respect  to  environmental 
perturbations,  as  well  as  angular  misalignment  of  the  fiber  probe.  Commercially  available  digital  signal 
processing  instrumentation  is  used  to  demodulate  the  interferometric  signal.  We  determined  that  the 

theoretical  dynamic  range  of  the  system  is  greater  than  10  . 

We  have  applied  this  interferometer  to  measure  the  transverse  deflection  of  a  MEMS  flexure  beam. 
We  investigated  the  response  of  the  flexure  beam  to  both  bipolar  and  unipolar  excitation  sinusoids  and  we 
characterized  the  tip  deflection  as  a  function  of  applied  voltage.  Our  experimental  results  were  in 
reasonably  good  agreement  with  coupled  BEA  and  FEA  simulation  of  the  flexure  beam. 

We  have  described  a  frequency-domain  algorithm  to  demodulate  an  optical  interferometer  using  DSP. 
We  have  shown  that  the  method  is  capable  of  unambiguously  detecting  vibrations  with  a  maximum 
amplitude  of  120  nm.  We  verified  this  frequency-domain  demodulation  technique  by  simulating  a  set  of 
interferometer  signals  and  by  characterizing  a  PZT.  While  time-domain  interferometric  demodulators  do 
exist,  our  frequency-domain  demodulation  scheme  manifests  certain  advantages  that  include  independence 
from  the  equilibrium  optical  path  length  difference,  the  reference  modulation  depth,  and  system  gain. 
Additionally,  because  the  target  frequency  must  be  greater  than  the  reference  frequency,  the  method  allows 
for  the  detection  of  very  high  frequency  vibrations,  limited  only  by  the  sampling  frequency  of  the  DSP. 

In  summary,  this  research  effort  has  produced  a  robust  interferometric  system  for  characterizing  the 
out-of-plane  deflection  of  MEMS  structures  and  a  DSP  algorithm  for  demodulation  of  the  interferometer. 
These  technological  advancements  yield  significant  progress  towards  the  development  of  smart  sensors 
based  on  interferometry  and  MEMS. 
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Appendix  A:  Quantum-noise-limited  sensitivity  of  an  interferometer 
using  a  phase  generated  carrier  demodulation  scheme 


At.  Introduction 

A  wide  variety  of  demodulation  techniques  have  been  developed  for  interferometric  sensors. 
Homodyne  demodulation  with  a  phase  generated  carrier  (PGC)  achieves  passive  operation,  large  dynamic 
range,  and  high  sensitivity.  Passive  operation  obviates  the  need  for  the  large  piezoelectric  phase 
modulators  and  fast  reset  circuitry  required  in  actively  stabilized  homodyne  systems  [Al],  In  addition,  a 

large  dynamic  range  on  the  order  of  109  has  been  realized  using  digital  signal  processing  to  perform  the 
demodulation  [A2].  In  this  paper,  we  derive  the  quantum-noise-limited  sensitivity  of  interferometers  using 
a  PGC  demodulation  scheme.  Quantum  noise  (also  called  shot  noise)  in  the  photodetection  process  sets  a 
fundamental  limit  in  the  performance  of  the  interferometric  sensor  [A3]. 

The  classical  Michelson  interferometer  is  the  optical  configuration  chosen  for  this  analysis.  However, 
the  following  sensitivity  analysis  is  applicable  to  all  two-beam  interferometers.  Figure  Al  shows  the 
system  layout.  We  form  the  PGC  by  dithering  the  position  of  the  reference  mirror  with  a  piezoelectric 
transducer.  This  interferometric  configuration  detects  changes  in  the  position  of  the  target  structure. 

A2.  Ultimate  System  Resolution 

The  PGC  is  imposed  on  the  interferometer  signal  by  sinusoidally  dithering  the  reference  mirror  at  an 
angular  frequency  of  coc  and  at  an  amplitude  of  -j,  where  A  is  the  light’s  wavelength.  The  reference 
beam’s  electric  field  can  be  defined  by 


Er  =  Ar  exp 


A 

jco0t-2k0Lr-2k0\  —  |sin(<yct) 


(Al) 


where  Ar  is  the  reference  field  amplitude,  coQ  is  the  light’s  angular  frequency,  ka  is  the  light  propagation 
factor  in  free  space,  and  Lr  is  the  distance  from  the  beam  splitter  to  the  reference  mirror.  The  signal 
beam’s  electric  field  can  be  defined  by 

Es  =  As  exp \jco0t  -  2 k0Ls  -  2 k08{t)\ ,  (A2) 

where  As  is  the  signal  field  amplitude,  Ls  is  the  distance  from  the  beam  splitter  to  the  target  structure,  and 
S(t )  is  the  surface  displacement  to  be  measured. 


Figure  Al.  Michelson  interferometer. 


Photodetector 
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If  we  assume  identical  polarization  states  of  the  two  fields,  then  coherent  interference  at  the 
photodetector  produces  a  resultant  electric  field  whose  magnitude  is  given  by 

^Resultant)  =  A'l  l  +  K~  +2K  C0s[2koJ(t)]  COS  [2kaA  -K  sin  a>ct] 

(A3) 

-  2 K sin[2k0<?(t)]sin[2k0A  -  ;rsin  coct\  }^2 

where  K  is  the  field  contrast  factor  ( As  l  Ar  )  and  A  is  the  quiescent  path-length  difference  between  the 
signal  and  reference  arms,  (Lr  —Ls  ).  To  express  the  limit  of  the  interferometer’s  sensitivity,  only  small 
displacements  are  considered  (i.e.,  kQS(t )  «  1 ).  Using  small  angle  approximations  yields 

|£Resultant|  =  4-[  1  +  K2  +  2K  cos(2/S;0A  -  ^sin  COct) 

(A4) 

-4Kk08(t)sin(2k0A-7Zsincoci)  J1'2 

The  last  term  within  the  square  brackets  of  Equation  (A4)  is  the  only  term  functionally  related  to  the 
surface  displacement,  8{t) .  This  is  the  only  component  of  the  resultant  field,  which  contributes  to  the 

signal  portion  of  the  photodetector  current,  is  .  The  entire  resultant  field  contributes  to  the  quantum  noise 
portion  of  the  photodetector  current,  i  . 

Now  that  we  have  defined  the  origins  of  signal  and  noise  in  the  system,  we  may  express  the 
performance  of  the  system  in  terms  of  the  SNR.  In  the  quantum  noise  limit,  the  photodetector  current’s 
SNR  is  given  by 


SNR  = 


(A5) 


where  (•)  indicates  a  time  averaging  operation.  With  respect  to  the  optical  detection  process,  the  signal 
current  is  defined  as 


=  P ' 


(A6) 


where  Ps  is  that  part  of  the  instantaneous  optical  power  incident  on  the  photodetector,  which  contributes  to 
the  signal,  77  is  the  detector’s  quantum  efficiency,  q  is  the  carrier’s  charge,  h  is  Planck’s  constant,  and  V 
is  the  light’s  optical  frequency.  Applying  the  signal  bearing  portion  of  Equation  (A4)  to  Equation  (A6) 
results  in 

is  =  C- ^  A Kk0S(t) sin(2.k0/S.~ 7t sin co ct)\,  (A7) 

where  C  is  a  constant  related  to  the  intrinsic  impedance  of  free  space  and  the  photodetector  area. 
Trigonometric  manipulation  and  use  of  the  Fourier-Bessel  expansion  reveals  the  mean-squared  signal 
current: 


=  C1At\J^- 
hv 


%K2klS2 


(0.-[l-J0(2^)cos(4koA)]. 


(A8) 
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To  determine  the  ultimate  system  resolution,  we  consider  the  quantum  noise  process  as  the  limiting 
noise  source.  The  detector’s  quantum  noise  is  related  to  the  time-averaged  photodetector  current,  (id  } ,  as 
follows  [A4] 


( 'ig)  =  2qAv{id ),  (A9) 

where  Av  is  the  detection  electronics  bandwidth.  This  limiting  mean-squared  noise  current  may  be 
expressed  as 


=  2qAv  ■  C ■  A2\  -  [  \  +  K2  +  2 KJ0 (n)  cos(2/L  A) 
h  v  ) 

-  4Kk08{t)J 0{n)  sin(2&0A) 


(A10) 


where  the  entire  resultant  field  of  Equation  (A4)  contributes  to  the  quantum  noise. 

The  SNR  is  found  by  substituting  Equations  (A8)  and  (A10)  into  Equation  (A5),  which  yields  the 
following  expression: 


SNR  = 


f  9  \ 

ca;?) 

1/2 

4K2  -4K2J0(2x)cos(4koA) 

hv-Av 

V  / 

1 +  K2  +  2KJ0{n)cos(2koA)  -  4 KkoS(t)J0(x)  sin(2/c0A) 

k08(t). 


(Ml) 


If  we  use  a  50/50  beam  splitter  and  assume  a  lossless  optical  system,  then  the  field  contrast  factor  is  unity 
and  the  SNR  becomes 


SNR  = 


(  m  f 

1  -  J0(2;r)  cos(4£0A) 

2  +  2J0(^)  cos(2&0A)  -4£oJ(t),/0(;r)sin(2A:oA) 

k08(t). 


(A12) 


where  is  the  laser’s  output  power.  This  is  the  quantum-noise-limited  SNR  of  an  interferometer,  which 
uses  a  PGC  demodulation  scheme. 

A3.  Results  and  Discussion 

Notice  that  the  quantum-noise-limited  SNR  is  a  function  of  A.  Figure  A2  illustrates  the  variation  of 
the  SNR  with  respect  to  the  quiescent  path-length  difference,  A ,  between  the  signal  and  reference  paths. 
Note  that  the  SNR  is  nonzero  for  practical  values  of  A .  This  result  is  in  contrast  to  homodyne 
interferometers,  which  do  not  use  a  PGC.  In  these  systems,  the  SNR  goes  to  zero  at  A  equal  to  an  integer 

multiple  of  y  [A3].  Active  stabilization  is  therefore  required  to  stabilize  the  system  at  quadrature  (i.e.,  A 

3 

equal  to  an  odd  multiple  of  y  ).  The  open-loop  interferometric  system,  however,  does  not  actively  stabilize 

the  system  at  quadrature.  Instead,  this  system  uses  a  PGC  to  extract  the  desired  amplitude  information  at 
any  path-length  difference,  A  . 

The  ultimate  resolution  of  the  system  is  calculated  directly,  using  the  results  of  the  SNR  analysis.  We 
set  the  requirement  that  the  SNR  equal  unity.  The  surface  displacement,  S(t) ,  then  assumes  its  value  at  the 

detection  limit,  which  we  denote  by  Amm  .  Let  us  assume  the  worst-case  scenario  for  the  mean  path  length 

3 

difference.  Figure  A2  indicates  that  minima  in  the  SNR  occur  at  A  equal  to  an  odd  multiple  of  y  .  If  the 
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Figure  A2.  Normalized  SNR  versus  path-length  difference 
between  the  signal  and  reference  arms  of  the 
interferometer. 


laser  power,  Pi ,  is  1  mW,  the  wavelength  of  light,  A  ,  is  0.6328  jtm,  the  detector  quantum  efficiency,  77,  is 

45%,  and  the  detection  electronics  bandwidth,  Av  ,  is  1  MHz,  then  the  minimum  detectable  displacement 
is  <Jmin  =  5  pm. 
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